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ABSTRACT 

The proposed global astrometry mission GAIA, recently recommended within the 
context of ESA's Horizon 2000 Plus long-term scientific program, appears capable of 
surveying the solar neighborhood within ~ 200 pc for the astrometric signatures of 
planets around stars down to the magnitude limit of V—17 mag, which includes late 
M dwarfs at 100 pc. 

Realistic end-to-end simulations of the GAIA global astrometric measurements 
have yielded first quantitative estimates of the sensitivity to planetary perturbations 
and of the ability to measure their orbital parameters. Single Jupiter-mass planets 
around normal solar-type stars appear detectable up to 150 pc (V ^12 mag) with 
probabilities ^ 50 per cent for orbital periods between ~2.5 and ~8 years, and their 
orbital parameters measured with better than 30 per cent accuracy to about 100 pc. 
Jupiter-like objects (same mass and period as our giant planet) are found with similar 
probabilities up to 100 pc. 

These first experiments indicate that the GAIA results would constitute an im- 
portant addition to those which will come from the other ongoing and planned planet- 
search programs. These data combined would provide a formidable testing ground on 
which to confront theories of planetary formation and evolution. 
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INTRODUCTION 



At t he I very end of 1995, ilia discovery ( Mayor & Qiielov 



and giant planets formed several AU from the ce ntral star, 
in contrast w ith the very short orbit al periods (|Mav or &j 



1995) df Che first Jupiter-mass (Mj) planet or biting a nui- 



mal star other than the Sun heralded the beginning of a new 
era of extraordinary discoveries in the realm of extra-solar 
planets, bringing with them the hope for a better under- 
standing of the formation and frequency of planetary sys- 
tems, and perhaps of bringing us closer to the ultimate goal 
of discovering extraterrestrial life. 

After four years since that discovery, spectroscopic pro- 
grams 1 
ets, i.e. 
cut-off 



Qucloz 1995) and high eccentricities ( Latham et al. 198!: 



Cochran et al. 1996; Mazch et al. 1996) found for several of 
the new discoveries. Their interpretation as bona-fide plan- 
ets rests on our understanding of correlations shown by their 
orbital and physical parameters, as recently disc ussed by 
Black (1997) and earlier by Duquennoy & Mayor (1991) in 
their work on solar-type binary stars. 

New models, which employ specific p hysical and dy- 



ave been able to reveal some tw e nty extra-solar plan- 
lqoo limit bolow tho 13 Mj 



objoota with a lowor i 

.vhich hac boon adoptod by Opponhoimor and Kulka 



namical mechanisms like in-sit u formation (Bodcnhcimer et 
al. 1999 ) or orbital migra tion ( |Lin et al. 1996 ; Trilling et al 



199S ; Murray et al. 199£ ) , have been proposed to justify the 



rni (199£) to differentiate giant planets from brown dwarfs. 

However, these discoveries have raised new and trou- 
bling questions in our understanding of the properties of 
planetary systems. The fundamental tenets upon which 
present theories are based include nearly circular orbits 
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presence of hot jupiters around normal stars, demonstrat- 
ing that the interplay between additional theoretical work 
and more observational data will be necessary for a contin- 
ued improvement in our theoretical understanding of how 
planets form and evolve, and where Earth-like planets could 
eventually be found. 

However, simply adding a few tens of additional detec- 
tions of giant extra-solar planets is not enough. A better un- 
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derstanding of the conditions under which planetary systems 
form and of their general properties requires large, complete 
samples of planets, with useful upper limits on Jupiter-mass 

planets at several AU from the central star. 

O ngoin g an d planned radial ve l ocity surveys ([Ma yor 



1996 



Queloz |1995 ;_ Cochran et al 

& Butljer 1992); [Marcy & Butler 1998p have started fill- 



Noycs et al. 1997| : |M arcji 



ing significant portions of the relevant parameter space. 
Searches based on relative astrometry from the ground and 
in space flGatcwuod 1987 ; Colavita ct al. 199£; Icolavita et 
al. 199^; |Mariotti et al. 1998j ; [Pravdo fc Shaklan 1996|) will 



be an important complement to the spectroscopic work and 
probably, the preferred means for establishing the existence 
of planets around young stars and that of low mass pla nets 



down t o a few Earth masses, as will be the case for SIM ( Bo- 
den et jtl. 1997|; [Unwin 1999|). 



A HIP PARC OS-Yike, space-borne global astrometric 
mission, which can survey the whole sky to faint magnitudes 
and with high astrometric accuracy, will enable the monitor- 
ing of large (> 10 ) samples of stars, with well understood 
completeness properties. This, depending on actual values 
of planetary frequencies ( Marcy et al. 1999| ), could yield the 
possibility of making firm measurements of statistical prop- 
erties of planetary systems. For, different correlations among 
orbital parameters (eccentricity, period or semi-major axis) 
and measurable differences in planetary frequency are likely 
to be generated by diverse planetary formation scenarios 
(core accretion and disk instability are the two known to 
date) and evolution mechanisms, as well as different forma- 
tion and evolution processes of the parent star (binarity, 
spectral type, metallicity, age) . An astrometric mission such 
as GAIA appears well poised for such a systematic census 
of planetary systems within ~ 200 pc from the Sun. 

The GAIA concept was originally proposed by Lin- 
degren and Perryman (1996) as a possible Cornerstone- 



class mission within the Horizon 2000+ program of scien- 
tific satellites of the European Space Agency. This satellite 
is designed to chart more than one billion objects (stars, 
extra-galactic objects, and solar system objects) on the sky 
down to the limiting magnitude of / = 20. The targeted fi- 
nal accuracy is ~ 10 fias on positions and parallaxes, and ~ 
10 /xas/year on proper mot ions at the reference magnitude 
of V = 15 for a G2V star ( |Gilmore et al. 1998| ), and for a 
mission life time of 5 years. 

In the following sections we show and discuss some rele- 
vant results derived from detailed end-to-end simulations of 
the data acquisition and analysis process for GAIA, which, 
as we will see, appears capable of discovering Jupiter-mass 



planets around r- 
earlier than K5). 



3x10 candidate stars (including dwarfs 



2 DATA SIMULATION 

The simulatio n cod e is an adaptation of that used by Gal- 
ligani et al. ( |l989| ) for the assessment of the astrometric 
accuracy of the sphere reconstruction in the HIPPARCOS 
mission . We generate catalogs of single stars randomly dis- 



tributec 1 on the sky; each run produces a sphere of 200 stars. 
As the satellite observing strategy (or scanning law) is most 
easily described in ecliptic coordinates, positions, proper 
motions and parallax (A, /3, fip, and n, respectively) 



are also given in the same coordinate system. For the mo- 
ment, /x/3 and 7r, as well as magnitudes and colors, are 
drawn from simple distributions which do not represent any 
specific Galaxy model; in particular, in each run the 200 
stars simulated have the same parallax, total proper mo- 
tions, magnitudes, and colors. 

The satellite is made to sweep the sky in such a way that 
the spin axis precesses around the Sun at a rate of about 6.5 
rev/year and with a constant angle of 43deg. Stars that at 
any given time are "seen" within a strip ~ 1 deg wide along 
the great circle (GC) being scanned are considered observed; 
a GC is completed in about 2.5 hours. Basic observations, 
which in principle can easily be derived from the fringe pat- 
tern measurements above, are the abscissae along a GC, as 
GAIA, like HIPPARCOS, makes very precise measurements 
in one dimension only. The mission lifetime is set to 5 years 
and the scanning law is such that the number of basic obser- 
vations per star is a function of ecliptic latitude: each given 
object on the sky is re-observed approximately every month, 
for a total of ~ 60 one-dimensional position measurements, 
on average. The minimum number of observations is ~ 30 
and occurs for stars at the ecliptic equator. The position of 
a star at a given time, as described by the combined effects 
of parallax and proper motions, is called here barycentric 
location and it has been described in Euclidean space (flat 
Lorentzian). General relativistic effects, which will have to 
be considered in the future (especially for the case of Earth- 
like planets), are not taken into account. 

Finally, gravitational perturbations (Keplerian mo- 
tions) induced by a (single) nearby orbiting mass are added 
to the barycentric location resulting in the "true" geometric 
location of a target. 

Simulated observations are generated by adding the ap- 
propriate astrometric noise to the true locations. The error 
sources considered in the simulation are discussed below. 

As in HIPPARCOS, GAIA will have two viewing direc- 
tions separated by a large angle named Base Angle (BA). 
It is by connecting directions far apart on the sky that the 
principle of space-borne global astrometry demonstrated by 
HIPPARCOS is implemented. Therefore, it is not surprising 
that the accuracy with which the BA is known through- 
out the mission is probably the most important single item 
within the GAIA concept. At the moment there are two 
mature optical designs for the BA. One configuration fea- 
ture two telescopes mounted on the same optical bench and 
pointing along the two different line-of-sights. The other 
design uses a beam combiner (an adaptation of that used 
on HIPPARCOS) which physically materializes the BA and 
multiplexes the two viewing directions into a single telescope 
unit; beam combiner and telescope are bolted to the same 
bench. The twin telescopes designed for the first configura- 
tion feature an off-axis monolithic primary (with no central 
obscuration), while the collecting telescope of the latter de- 
sign is a Fizeau interferometer with a baseline of 2.45m. 

The details of th e GAIA optical configuration (see e.g. 
Gilmore et al. (1998), and references therein) are not crit- 



ical to our discussion of photon-driven astro metric errors. 
The ma jor difference with our earlier work (Casertano et 



al. 1995) is that both designs feature a significantly increased 



telescope aperture. The monolithic configuration has a rect- 
angular primary of 1.7m by 0.7m, and each circular aperture 
of the interferometric option is 0.65m in diameter. High- 
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Table 1. Photon and total error for a single observation. This 
is the combination of ~ 8 consecutive field- of- view crossings, for 
any given direction on the sky, each comprising 20 elementary 1- 
sec exposures, for a total duration for the single integration of ~ 
160 sec. These values strictly apply to the inter ferometric option 
discussed in the text and to a G2V star and negligible interstellar 
absorption (Av = U). The pure photon noise values increase by 
~ 30 % fwi an AOV siwi — and imp'iuve by i*= 20 % fwi a K3V. 
More explanations are provided in Sections 2.2-2.3 



V mag 


Photon error 


Total error 


( ^tas) 




( fias) 


(<r b = 30 pm) 


(°b = 10 pm) 


10 


2.4 


10.2 


3.8 


12 


5.9 


11.5 


6.6 


15 


24.6 


26.5 


24.8 


17 


70.3 


71.0 


70.4 


18 


129.9 


130.3 


129.9 



accuracy measurements of the position (phase) of the PSFs 
are made directly on the focal plane using CCD detectors 
(see next section). 



2.1 Payload configurations 

The payload design has greatly a dvan ced since the idea 
sketched in Lindegren & Per ryma n (1996) and further devel- 
oped in Loiseau & Shaklan (199C). The two industrial stud- 
ies commissioned by ESA have looked into t he feasibility of 
two different options for the GAIA payload ( |Gilmore et al 



199S 



an all-passive configuration with two identical tele- 
scopes (with rectangular-shape monolithic primary mirrors) 
to be operated in L2, and an all-active configuration with 
an interferometric (diluted) beam combiner and a Fizeau 
interferometer as light collector behind it to be operated in 
geosynchronous orbit. Stability of the BA is passively main- 
tained in the case of the monolithic configuration by uti- 
lizing a sophisticated active thermal control system (which 
must operate at the fiK level) and a silicon carbide structure 
for the optical bench. On the other hand, the interferomet- 
ric configuration is all-active, with the stability achieved by 
real-time monitoring implem ented via high-accuracy laser 
metrology (Gai et al. 1997) and control of all critical de- 
grees of freedom. 

As we are more familiar with the interferometric op- 
tion for GAIA (Lattanzi et al. 1997), we will be referring to 
that in the discussion below. However, precision and accu- 
racy requirements and estimates are quite similar for both 
configurations, especially at the bright magnitudes we are 
concerned with in this paper. 



2.2 Photon noise and point spread function 
measurement 

The ability to measure accurately the position of a star 
ultimately depends on the width and shape of the point 
spread function (PSF) of the imaging system and on the 
number TV of photons detected. F or a well-sample d PSF, 
the theoretical limit is shown by (Lindegren 1978) to be 
e = \/(4tvx t 



N), where is the rms size of the aper- 
ture in the measurement direction. For two circular aper- 
tures of diame ter D and with a central separation B, we 
have a; rms = y/(B/2) 2 + (D/4) 2 ; for the baseline GAIA pa- 
rameters (73 = 2.45 m, D = 0.65 m), x rms = 1.24 m (Lin- 



degren fc Perryman 1995] ), For A = 750 nm (the baselined 
operational wavelength), this translates into a theoretical 
monochromatic measurement accuracy of lOmas/x/TV- 

The measurement accuracy is degraded for non- 
monochromatic measurements, by about a factor 2 for a 



Gaussian filter centered at 7500 A and 2000 A wide (Gai et 



al. 1998). In addition, the requirement of optimal sampling 
may be difficult to achieve, since the central fringe is only 
about 0'.'04 wide and the field of view is 1 deg. In practice, 
this will probably cause a loss of accura cy of about 20-40 
per cent ( |Gai et al. 1995 ; Gai et al. 1998 ). In the following, 
we assume a "best reasonable" single-measurement photon 
noise error of 24mas/\/7V : . 

Since scans overlap partially, each "observation" of a 
star will consist of about 8 consecutive scans with 20 sec of 
integration time allocated per scan, a total integration time 
of 160 s. Assuming a total collecting area of 0.664 m 2 (2 
apertures of 0.65 m diameter) and a total system efficiency 
of 20 per cent, a star with V = 15 would generate about 10 6 
photons per observation, corresponding to a photon-limited 
measurement accuracy of the fringe position of 24 /j,a.s. The 
accuracy scales with the inverse square root of the flux. This 
approximate calculation agrees with the values in Table |l[ 
obtained for the current version of the interferometric op- 
tion. The photon noise values (second column in Table Q) for 
10 ^ V ^ 18 stars were obtained with 3.1 pixels per fringe 
period, a RON of 3 electrons/pixels, and a DQE of 0.6. Ob- 
jects brighter than V = 10 mag will be also observed by 
GAIA. However, limitations on the dynamic range achiev- 
able with CCDs, requires to limit the exposure time for the 
brighter stars. Thus is practice, the precision starts to level 
off brighter than V ~ 10. 

This accuracy is based on the photon statistics only, 
and does not take into account possible systematic effects, 
such as distortions in the optical system or in the detector, 
imperfection in the fringes, and so on. Many such systematic 
effects can be calibrated using closure methods, others will 
require enhancements in system design. 



2.3 Accuracy of the Base Angle 

The systematic effects mentioned above will lead to resid- 
ual (systematic) errors in the knowledge of the BA which 
do not scale with magnitude. Ultimately, these will be the 
errors which will limit the maximum accuracy of GAIA for 
bright sources. We combine such residual errors in what is 
hereafter called residual bias, or simply bias. In the present 
error budget we assume that the bias can be described as a 
stationary stochastic process; therefore the bias contributed 
at the single-observation level to each object (Table ^) scales, 
like photon noise precision, with the average number of ob- 
servations collected over the mission lifetime. 

The BA can be measured and monitored accurately over 
time intervals longer than a full revolution (^ 3 hours) by 
making use of the 2-7T closure properties of the consecutive 
scans. This is an especially important feature of GAIA whose 
potential, at the mas level, has been beautifully proven by 
HIPPARCOS. However, the bias over shorter time scales 
need to be controlled by ensuring that the relative positions 
and shapes of all optical elements of the beam combiner 
do not vary throughout the observation. The necessary sub- 
nm accuracy will probably be achieved by a combination 
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of passive control and of laser metrology. For simplicity, we 
summarize the overall accuracy (bias) with which the (rel- 
ative) positions of the mirrors of the beam combiner must 
be known by a single number, the "baseline error" <7b- This 
helps visualize the complex interplay of the beam combiner 
mirrors (which materialize the two baselines of the interfer- 
ometric design) with a more familiar bias of the type B x S, 

where B is the baseline and 5 the angular uncertaintyPj. This 
also helps understand the bias contribution in Table Til For 
example, if the photon error is subtracted (in quadrature) 
from the total error for the case with Ob = 30 pm, the re- 
sulting value, ~ 10 /Lias, represents the angular bias cor- 
responding to that linear "baseline error". Therefore, the 
requirement on the metrology is ~ 3.8 times more stringent 
that one would have guessed from the simple calculation 
B x S, with B=2.45 m and S=10 /ias. This example shows 
that the baseline bias is driven by the dimension of the single 
mirrors (D=0.65 m) forming the baseline (B/D~3.8) within 
the beam combination system. 

Laser metrology in laboratory settings has already 
achieved extremely high perfor mances, with relative mea- 



sure ments at the picqmeter level ( Siirsel 199E ; Noecker et al 
1993; N oecker 1995 ; Reasenberg et al. 1995 ). Such precision 
has been reached over short (few wavelengths) variations in 
path length, which are appropriate to the GAIA design if a 
good active thermal control is included. 

However, the few pm error quoted refers to the precision 
and stability of a one-dimensional laser gauge measurement 
of a single optical path. Maintaining the accuracy of the in- 
terferometer baseline is much more complex, first, because 
the three-dimensional positions of several optical elements 
may need to be monitored simultaneously, and second, be- 
cause of the possible differences between the optical path 
of starlight and of the laser gauge beams. Noecker ( 1995 ) 
lists a number of possible systematic errors for the POINTS 
mission concept. An experiment to demonstrate picometer 
laser metrology for a three-dimensi onal system on the GAIA 
scale is underway (Sai et al. 1997). For the moment, it re- 



mains difficult to give firm figures for the baseline accuracy 
that will eventually be achieved. In Table jj], we consider two 
cases which probably bracket realistic expectations: a more 
"conservative" accuracy o"b = 30 pm, which already appears 
within reach from the preliminary results available, and an 
"optimistic" accuracy <Tb = 10 pm. Notice that these bias 
values are intended at the level of what is called here the 
single-observation error. The accuracy levels at the end of 
the mission will be ~ 3 pas and 1 j^as for Ob = 30 pm and 
10 pm, respectively, if errors in successive passes are uncorre- 
cted, as discussed previously. Notice that, for Ob = 30 pm, 
the astrometric precision begins to level off brighter than 
V = 12. 

Finally, it is important to mention that the numbers 
quoted for the linear bias apply to a 3-dimensional monitor- 
ing of the relevant structure. The laser beams will have to 
monitor the corresponding one-dimensional degrees of free- 



t While the main contribution to the baseline error will probably 
come from the relative position of the beam combiner mirrors, the 
motions and/or distortions of other optical elements (such as the 
primary apertures of the Fizeau interferometer) can also generate 
an effective baseline error. 



dom with significantly better resolution — a factor < 2 for 
the current design. 



3 DETECTION AND ORBIT 

DETERMINATION METHODS 

The magnitude of the gravitational perturbations induced 
by a planet on the parent star, as seen by an astrometric 
mission, can be quantified through its astrometric signature 
a defined as: 



Mp a P 
M s d 



(1) 



where M p , M s are the masses of the planet and star re- 
spectively, o p the semi-major axis of the planetary orbit, d 
the distance of the system from us. If a p is in AU and d in 
parsec, then a is expressed in arcsec. 

GAIA's sensitivity to this signal depends of course on 
the errors of each measurement, with theoretical values 
listed in Table |lj We have verified through tests with differ- 
ent values of (see section 4.1) that, as could be expected, 
the detection probability depends on and a only through 
their ratio, the "astrometric" signal-to-noise ratio 



S/N = a/ cry, , 



(2) 



so that the results obtained can be straightforwardly 
rescaled to different measurement errors. 

For simplicity, we have thus kept the single-observation 
error fixed at — 10/ias throughout our simulations. This 
is the value expected for relatively bright stars (V ^ 12 
mag, corresponding to the Sun at 200 pc), with a conserva- 
tive baseline error of 30 pm (see Table 111) . Our simulations 
should thus give realistic estimates of: a) GAIA's detectabil- 
ity horizon of planetary mass companions to solar-type stars 
in the vicinity of our Sun, and b) limits on distance for ac- 
curate orbital parameters determination. 

The starting point for our two-level statistical investi- 
gation is the computation of detection probabilities, which 
will in principle depend upon 1) mission parameters, 2) noise 
sources, and 3) orbital elements. The contributions to points 
1) and 2) are listed in Section 2. As for point 3), we will ex- 
press the detection probabilities as function of the period 
P and the signature a, which we expect to be the major 
contributors, and generally average over the expected dis- 
tribution of the other orbital parameters. 



3.1 Detection 

Our first detection method applies a standard \ 2 test (with 
the confidence level set to 95 per cent) to the residuals ip—ip r , 
where the xj) are the actual measurements, and the ijj r are 
the GC abscissae recomputed after fitting the observations 
of each star with a single-star model, i.e., solving only for 
the five astrometric parameters appropriate for a single star 
(position, parallax, and proper motion). Being Xo the value 
provided by the null model (no planet), the test fails when 
Pr(x 2 < Xo) > 0.95. In this case the planet is considered 
detected. Note that this method only measures deviations 
from the single-star model, it makes no assumptions on the 
nature of the residuals nor does it provide initial guesses for 
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the computation of the planet's mass and orbital parame- 
ters. 

We employed the x 2 test to analyze 160 000 stars uni- 
formly distributed on the sky, perturbed by dark compan- 
ions inducing astrometric signatures a ranging from 10 to 
100 fias, with periods P between 0.5 and 20 years. The re- 
maining orbital elements were distributed randomly in the 
ranges: 0° < i < 90°, < e < 1, < Q. < 2tt, < u < 2tt, 

^ r ^ P. The single observation error was set to = 
10 /ias, thus the signal-to-noise ratio varied in the range 

1 ^ S/N ^ 10. Finally, we repeated the same simulations 
without planets, to verify the correct behaviour of the test 
and the choice of the confidence level. As expected the num- 
ber of false detections was ~ 5 per cent. 

3.2 Orbital parameters 

Once a planet is "detected", the goal is to derive reliable 
estimates for of its orbital elements and mass. For a complete 
reconstruction of the orbital geometry, we implemented an 
analytic model in which the observation residuals ^) — %p r are 
evaluated with a recomputed abscissa of the form: 

ipr = 4>r{\ P, H\, fip, tv, Ti, T 2 , T s , T 4 , e, P, r) (3) 

The 1-d measurement along the scanned GC is thus 
expressed as function of both the standard astrometric pa- 
rameters (A, f3, n\, fip, 7r) and the parameters describing 
the stellar relative orbit with respect to the barycenter of 
the planetary system: period P, periastron epoch r, eccen- 
tricity e and the 4 Thiele-Innes elements T, functions of 
semi-major axis a, inclination i, argument of periastron ui, 
position angle of the ascending node Q. 

The solution of the non linear system of observation 
equations is obtained employing an iterative linear Least 
Squares procedure by means of which the entire set of orbital 
elements can be simultaneously estimated within a well de- 
fined accuracy level. The details of this procedure are given 
in the following sections. 



4 RESULTS 

In this section we account for what is our present under- 
standing of GAIA's sensitivity to astrometric perturbations 
induced by Jupiter-mass planets orbiting around nearby 
stars. Given the present evidence of a number of extra-solar 
planets, we also present results on how GAIA would perform 
on a selection of three such systems. 

4.1 Detection probabilities 

Figure ^ gives the percentage of failure of the \ 2 test on the 
single-star hypothesis in the case of a^, = 10 fias. We note 
that at relatively low S/N ratios the detection probability is 
dominated by sampling of the orbital period. At higher S/N 
values orbital sampling is less critical and long period plan- 
ets (up to about twice the mission duration) are detectable: 
as a matter of fact, when S/N — > 10 detection probability 
reaches about 100 per cent. 

The dip at P = 1 year was somewhat expected, as both 
orbital motion and parallactic factor have the same period. 
However, the decrease is small, indicating that the coupling 




Figure 1. Planet detection probability (percentage of failures of 
the x 2 test) as function of the astrometric signature a and of the 
orbital period P, for =10 fias. The percentage of detection of 
each point is based on 200 random planetary systems uniformly 
distributed on the sky. 
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Figure 2. Iso-probability contours (solid lines) for 25 per cent, 
50 per cent and 95 per cent of detection probability, compared 
with Kepler's third laws (dotted/dashed lines) for systems with 
Jupiter-Sun masses at D=50, 75, 100, 150 and 200 pc. Jupiter- 
like planets (P = 11.8 yr) appear detectable, with probability ^ 
50 per cent, up to a distance of 100 pc (vertical line). 



is less critical than might have been anticipated. This is 
probably due to the fact that the parallactic motion has 
fixed phase and aspect ratio, and a relatively small mis- 
match in any of the orbital parameters — phase, inclination, 
eccentricity — is sufficient to separate the two signals. 

Figure |l] also shows that the x 2 test quickly loses its 
sensitivity for S/N approaching unity; thus planets for which 
the error in individual observations is comparable to the 
apparent semi-major axis are essentially undetectable with 
this technique. 

As mentioned before, the results of Figure ^ apply to 
measurement errors other than the assumed = 10 />as, 
as long as the S/N value remains the same. For example, 
for a single observation error = 1 /^as, the detection 



© 2000 RAS, MNRAS 000, 0-|l| 



6 M.G. Lattanzi et al. 

probability is exactly the same as shown in Figure |l| but for 
amplitudes ten time smaller. 

We can get important physical informations from the 
statistical results of Figure nj, simply changing the interpre- 
tative perspective from which we are leading the discussion. 
The solid lines in Figure Q are the empirical relations, de- 
rived from Figure [l], that express the amplitude of the per- 
turbation as a function of orbital period needed for detec- 
tion probabilities of 25, 50, and 95 per cent, respectively. If 
the orbital period is shorter than the assumed mission life- 
time of 5 years, the detection probability is about 50 per 
cent for S/N ~ 1. For longer periods, the sampling of the 
orbit is worse, probability drops significantly and a much 
higher signal is required for the planet's signature to be de- 
tected. T his is in qualitative agreement with the results of 
Babcock (1994), who studied the detection and convergence 
probability of a complete orbital model for simulated plane- 
tary systems as observed by the mission POINTS. Babcock 
did find a slightly larger sensitivity on planet's period, man- 
ifested in an earlier turn-up and steeper slope at long pe- 
riods of the 50 per cent probability curve; this most likely 
depends on the fact that the determination of reliable orbital 
elements is more challenging than detection only. 

Let us now discuss the other elements of Figure [| The 
overplotted dashed and dotted lines represent the signal ex- 
pected for a Jupiter-mass planet and solar-mass star at var- 
ious distances and orbital periods, obtained by substituting 
Kepler's third law in the expression for a: 



O.OSeSl.O 



O.OSeSl.O 



M K P 2/3 



M, 



2/3 d 



(4) 



The vertical solid line at P = 11.8 years indicates 
the locus of the actual Jupiter-Sun system at different dis- 
tances. Jupiter-Sun systems appear detectable with proba- 
bility 50 per cent up to a distance of 100 pc, while Jupiter- 
mass planets with shorter periods can be detected to larger 
distances: to 150 pc for periods between 2.5 and 8 years. The 
relation can be rescaled to lower planet masses by reducing 
the distance of the system in proportion to the planet mass, 
since it is the ratio M p /d that enters in the astrometric sig- 
nature a. 

Periodic photospheric activity across the disk of the par- 
ent star (e.g. star-spots cycles) can induce displacements in 
the position of the photocenter, thus adding astrometric sig- 
nal which could be difficult to disentangle from the planet 
signature. However, the magnitude of such intrin sic astro- 
metric noise appears to be at most of a few /ias (Woolf & 



Angel 1999) for a solar-type star at a distance of 10 pc, i.e., 
significantly smaller than the expected single-measurement 
errors. 

Finally, although the single-observation error signifi- 
cantly deteriorates with magnitude (Table , jovian planets 
around late type stars can reliably be detected at relatively 
large distances. 



4.2 Orbital parameters estimation 

The next step beyond the simple detection of candidate 
planets is the task of estimating the orbital parameters of 
each system and the mass of the unseen companion. 

Our method, described in more detail below, consists 
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Figure 3. Convergence probability to a 10 per cent accuracy for 
a, P, e, i, as function of the distance from the Sun. 



of applying an iterative non-linear least-squares fitting pro- 
cedure to each of the simulated orbits; the 'known' orbital 
parameters are utilized as initial guesses to start the fit. 
Convergence of the non-linear fitting method and quality of 
the orbital solutions can both be significantly affected by the 
choice of the starting guesses. Therefore, the use of the true 
values of the orbital parameters to initialize the fit leaves 
open some important questions related to how and to what 
extent effective starting values, i.e., leading to successful or- 
bital solutions, will be identified from the data as function of 
actual performances of the satellite, uncertainties in the er- 
ror model, and properties intrinsic to planetary systems. For 
this, realistic global search strategies must be implemented 
and double-blind test campaigns conducted. Work on these 
issues is in progress and will be reported in the future. 
Instead, this work focuses on the important goal of gaug- 
ing GAIA's ultimate ability in measuring planets (properly, 
single giant planets orbiting single solar- type stars). An effi- 
cient way of achieving this is by assuming perfect knowledge 
of GAIA's characteristics (mainly, error model and satellite 
attitude) and, indeed, by using the known values of each 
orbital parameter as suitable initial guesses for the least- 
squares solutions Then, the post-fit differences to the true 
values of the parameters should be a reliable measure of 



t Note that, as explained later in this section, the fitting program 
does not know that the initial values provided for the parameters 
are also their true values 
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Figure 4. Same as Figure for convergence to a 30 per cent 
accuracy level. 
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Figure 5. Same as Figure ^, but for convergence to a 50 per cent 
accuracy level. 



the accuracy in orbit reconstruction that can ultimately be 
achieved by GAIA with the given measurement errors. 

Similarly to what we did in section [O], we look at the 
accuracy of the results primarily as function of distance and 
orbital period. We consider systems with Sun- Jupiter masses 
and three values of the period: 0.5 years, to test GAIA's 
ability to cope with poorly-sampled motion; 5 years, a near- 
ideal case where the orbital period is as long as the mission; 
and 11.8 years, the true period of the Jupiter-Sun system, 
which stretches GAIA's ability to solve long-period orbits. 
For each case, we generate 200 systems randomly placed on 
the sky, with randomly distributed values of the other orbital 
parameters. These 200 systems are then placed all at the 
same distance, and the simulation is repeated for distance 
values ranging from from 10 to 100 pc. 

A fit to the simulated observations is made directly at 
the GC level, taking into consideration a theoretical model in 
which GAIA's unidimensional measurements are expressed 
as functions of the 5 astrometric parameters plus the 7 
orbital elements of the keplerian orbit: for the latter we 
adopted a slightly linearized analytic form where semi-major 
axis, inclination, periastron longitude and position angle of 
the ascending node are combined in the four Thiele-Innes 
elements, while only orbital period, eccentricity and perias- 
tron epoch are left unchanged. 

The first solution of the linearized system is thus found 
with respect to the initial guess. Then, the linearization at 
step k is updated with respect to the solution obtained at 
step fe— 1 and the process repeated until the differential cor- 



rections Sa^k to each of the 12 parameters satisfy the rela- 
tion: \5ai,k/cti,k-i\ < 10~ 6 , where i = 1,2, . . . , 12 and a^k-i 
are the parameters adjusted at the previous step. We then 
record what fraction of the final values for each parameter 
falls within a certain fractional error of the true value: this 
we call the "convergence probability", which is a function 
of the distance of the system, the period, and the desired 
fractional error. 

We evaluate the convergence probability for the param- 
eters which are most likely to affect the efficiency of the 
reconstruction of the orbit and of the mass determination, 
namely semimajor axis a, period P, inclination i, and eccen- 
tricity e. For each parameter, we consider fractional errors 
of 10, 30, and 50 per cent. 

The results are shown in Figures ^ to |j| Points of dif- 
ferent shapes correspond to different periods: triangles for 
0.5 years, diamonds for 5 years, and crosses for 11.8 years. 
Figure ^| shows the probability of convergence within 10 per 
cent of the true values; Figure ^ to within 30 per cent; and 
Figure ^ to within 50 per cent. In each Figure, different pan- 
els correspond to different quantities. 

As expected, the 5 year-period case is the best of the 
three. The 11.8 year period, although it corresponds to a 
larger signal (as the semimajor axis, and thus a, increases 
with period), suffers from the incomplete sampling of the 
orbit during the mission life-time, while the short-period 
case suffers from both the smaller signal amplitude and a 
(generally) non-optimal timing of the observations. 

The periodicity of the signal is the characteristic which 
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Figure 6. The effect of eccentricity: the plots show the 20 per 
cent accuracy level in cases of Sun- Jupiter systems with orbital 
period P = 0.5 years and e set at: 0.1 (diamonds), 0.3 (triangles) , 
0.5 (stars), 0.7 (points), 0.9 (squares). 



can be evaluated with the best accuracy: P is the only el- 
ement which is always estimated with high accuracy (bet- 
ter than 10 per cent) throughout the ranges covered by our 
simulations. We also note that the semimajor axis a and the 
inclination i behave similarly to each other. This may be in 
part a consequence of the use of the Thiele-Innes representa- 
tion, since both parameters are obtained from combination 
of the Thiele-Innes elements. 

Somewhat surprisingly, orbital eccentricity — which in 
the previous simulation is assumed to be distributed uni- 
formly from to 1 — has a very significant effect on the qual- 
ity of the estimated orbital parameters, including the mass of 
the planet. We illustrate this in Figures ^| through^, where 
the different symbols now refer to different orbital eccen- 
tricities, and the convergence probability is given for a fixed 
fractional error of 20 per cent. Different figures correspond 
to different orbital periods. Generally, high eccentricities de- 
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Figure 7. Same as Figure for Sun- Jupiter systems with P = 
5 years. 

teriorate convergence percentages for all orbital parameters. 
The reason is that the regularly-spaced observations of a 
survey satellite, such as GAIA, may be ill-suited to sample 
an orbital motion with large velocity variations, as happens 
for high eccentricities. The effect is especially prominent for 
the long-period case, P — 11.8 years, for such orbits will 
often never be observed during the periastron. 

Our main results can thus be summarized as follows: 

- For a uniform eccentricity distribution in the range 
^ e ^ 1, then more than a half of all existing Jupiter- 
mass planets orbiting around solar-type stars with periods 
comparable to the mission lifetime can be detected, and their 
masses and orbital parameters can be estimated to an accu- 
racy of 10 per cent up to distances of about 100 pc. From the 
completeness limit of the HIPPARCOS catalog (V ~ 7.5-8 
mag) and considering spectral types no later than G5, we es- 
timate 20 000, 65 000, and 150 000 stars to 100, 150, and 200 
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Figure 8. Same as Figure [ij, for Sun- Jupiter systems with P 
11.8 years. 



pc respectively. In this estimate the contribution from early 
type stars and giants is negligible. Note that these numbers 
increase significantly when considering spectral types earlier 
than K5: calculations based on current Galaxy mo dels pre- 
dict > 300 000 stars within 200 pc from the Sun (Lattanzi 
et al 



11)99) 



- For a fractional error no worse than 20 per cent, the 
distance limit is reduced to 50 pc for short-period and long- 
period orbits (0.5 and 11.8 years); 

- Low-eccentricity systems are easier to detect and solve: 
this could indicate a potential limit in GAIA's capability of 
reconstructing the precise behaviour of the mass function 
of low mass companions to solar-type stars in the neigh- 
borhood of our Sun, as to our knowledge first observational 
data seem to indicate than brown dwarfs are more likely 
to revo lve around t he parent stars on significantly eccentric 
orbits (Black 1997). On the other hand, low orbital eccen- 



tricities may be prevalent for longer periods, thus enhancing 
GAIA's capabilities for such systems. 
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Figure 9. Comparison between GAIA 50 per cent probability 
curve (solid line) and the minimum astrometric signature (a) 
induced on the parent star by presently known planetary mass 
objects (0.5 Mj S^. M ^ 11 Mj ) in the neighborhood of our sun, as 
function of orbital period (log-log plot). The dahsed, dotted, and 
dashed-dotted lines represent the astrometric signatures induced 
on a solar-like star by the major solar system bodies, at increasing 
distance from the observer 



4.3 47 Uma, 70 Vir, and 51 Peg 

We now consider more specifically how the known plane- 
tary systems fall within GAIA's capabilities. Figure ^ shows 
period and astrometric signature for the planets detected to 
date by radial velocity measurements as compared to GAIA 's 
50 per cent iso-probability curve. This is the minimum sig- 
nature, corresponding to an orbital inclination of 90°; the 
radial velocity method cannot generally remove the degen- 
eracy between planetary mass and inclination. Although all 
systems are very close to the Sun (d < 40 pc), several plan- 
ets have very small signatures due to their short orbital 
period (and thus small semi-major axis, by Kepler's third 
law). Some of the shortest-period objects, i.e, the 51 Peg- 
class planets, may prove difficult to detect and measure with 
GAIA. On the other hand, planets with signature ^ 100 /ias 
will be easily detected, and their orbital elements can be 
found with high accuracy by GAIA. 

The considerations and results obtained in the previous 
sections from a statistical viewpoint are obviously prelimi- 
nary to the much more difficult task of the development of 
self-consistent detection and orbital parameters estimation 
algorithms devoted to the investigation of individual objects, 
in which a more complete and realistic error model can be 
taken into account: to this end detailed system studies are 
in progress and will be presented elsewhere. It is neverthe- 
less of interest a first glance to how well GAIA will behave 
once at work: thus, to test further GAIA's capabilities of de- 
tecting periodic oscillations and signatures due to planetary 
companions around stars in the neighborhood of our solar 
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Figure 10. Planetary perturbations of 47 UMa as "seen" by GAIA for an orbital inclination of 45° ■ The top left panel shows the path 
(solid line) of the star over a 5-year period. The * represent the true location of the star at the time of the simulated GAIA observations. 
The upper right panel zooms on the 0.1" X 0.1" region marked with a small square in the first panel. As for the two bottom graphs, 
they each represent a different 1-mas 2 enlargement of the top right panel. Again, the solid line represents the true (simulated) path of 
the parent star of the system. A single-star fit for 47 UMa produces the dotted path with crosses marking the calculated positions at 
the same observing times as before. The dashed line shows the recomputed positions after fitting the observations with a 10-parameter 
model (five astrometric parameters plus the four Thiele-Innes elements and the orbital period) 



system, we have concentrated our attention on the results 
obtained with radial velocity techniques, for what concerns 
three of the star-planet systems known to date: 47 UMa, 70 
Vir and 51 Peg. All these stars are nearby and very similar 
to our Sun. 



Accordi ng to the spectroscopic measurements ( Butler & 



Marcy |L996| ; |Marcy fc Butler 1996| ; [Mayor fc Queloz 1995[ ) 



around these three stars orbit planetary bodies of minimum 
mass, respectively, ~2.46 Mj, ~6.50 Mj, and ~0.5 Afj, with 
orbital periods of about 3 years, 4 months, 4 days. This 
translates in the following astrometric signatures: 

a (47 UMa) > 361.88 /xas 

a (70 Vir) > 167.93 fias 

a (51 Peg) 1.66 fias 

Because of its extremely short period, the astrometric 



perturbation in the case of the 51 Peg system is so small 
that reliable estimates of its astrometric orbit with GAIA 
will be highly difficult. On the other hand, the detection 
of the signatures induced on 47 UMa and 70 Vir will be 
extremely easy for a GAIA-like satellite. 

We have generated 100 simulated planetary systems on 
the celestial sphere, respectively identical to 47 UMa and 70 
Vir, assuming perfectly circular orbits and stellar masses M 
= Mq. The unknown orbital inclination was chosen to be 
i = 45°. 

First of all, we applied the \ 2 test to the simulated 
observational data, with a conservative single-observation 
error o$ = 10 ^tas. The detection probabilities derived are, 
as expected, always about 100 per cent, given the signal-to- 
noise ratio to be always in the regime S/N 3> 1. 

In order to obtain a visual indication of the quality of 
the reconstruction, we then compared the apparent motion 
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Table 2. Results of fitting GAIA observations of 47 UMa: the true orbital period simulated is Ptrue = 2.98 yr, 
the star's parallax is irtmc = 71.03 mas. The estimated values are averages over 200 simulations per orbital 
inclination bin. The RMS errors are derived from direct comparison to the true values. 



i 


atruc (£tas) 


Cfit (/ias) 


Pflt (years) 


7T flt (mas) 


M p ,truc (M.j) 


Mp, flt (Mj) 


20° 


1058.08 


1060.92 ± 15.15 


2.951 ±0.125 


71.0377 ±0.0026 


7.19 


7.27 ±0.24 


30° 


723.77 


725.56 ± 14.45 


2.956 ±0.142 


71.0377 ±0.0029 


4.92 


4.97 ±0.19 


40° 


562.99 


565.62 ± 14.97 


2.955 ±0.119 


71.0381 ±0.0027 


3.83 


3.87 ±0.15 


50° 


472.41 


470.75 ± 14.82 


2.958 ±0.158 


71.0377 ±0.0030 


3.21 


3.22 ±0.15 


60° 


417.87 


419.63 ± 17.07 


2.959 ±0.179 


71.0372 ±0.0028 


2.84 


2.87 ±0.16 


70° 


385.11 


384.36 ± 14.48 


2.973 ±0.112 


71.0378 ±0.0026 


2.62 


2.62 ±0.12 


80° 


367.47 


367.93 ± 16.38 


2.962 ±0.157 


71.0381 ±0.0025 


2.50 


o r. 1 _l_ n i a 


90° 


361.88 


362.64 ± 15.17 


2.930 ±0.159 


71.0377 ±0.0023 


2.46 


2.50 ± 0.14 


Table 3. Results 


of fitting GAIA observations of 70 Vir: the true orbital period simulated 


is Ptruc = 0.32 


the star's parallax 


is 7rtruc = 55.22 mas. 








i 


Otrue (/^as) 


ffit (A*as) 


Pflt (years) 


7T flt (mas) 


M Pl true (Mj) 


M p ,flt (Mj) 


10° 


967.05 


974.56 ± 15.83 


0.3198 ±0.0015 


55.2216 ± 0.0033 


37.43 


37.73 ± 0.62 


20° 


490.98 


495.70 ± 15.54 


0.3199 ±0.0016 


55.2211 ±0.0027 


19.00 


19.19 ±0.60 


30° 


335.85 


339.64 ± 14.36 


0.3201 ± 0.0021 


55.2208 ± 0.0030 


13.00 


13.14 ±0.56 


40° 


261.25 


265.61 ± 15.37 


0.3201 ± 0.0023 


55.2215 ± 0.0030 


10.11 


10.28 ±0.60 


50° 


219.21 


218.48 ± 15.39 


0.3198 ± 0.0023 


55.2209 ± 0.0030 


8.49 


8.46 ± 0.60 


60° 


193.91 


195.89 ± 17.03 


0.3199 ± 0.0029 


55.2205 ± 0.0028 


7.51 


7.58 ±0.66 


70° 


178.70 


181.71 ± 14.71 


0.3198 ±0.0032 


55.2215 ± 0.0025 


6.92 


7.03 ±0.57 


80° 


170.52 


172.00 ± 14.23 


0.3197 ± 0.0030 


55.2211 ± 0.0026 


6.60 


6.66 ±0.55 


90° 


167.93 


169.69 ± 14.31 


0.3198 ± 0.0042 


55.2210 ± 0.0029 


6.50 


6.57 ±0.56 



of the star with the motion obtained with a single star fit 
and that obtained by a fit with the fitted parameters for the 
planet. Figure ^ shows the differences in 47 UMa's path on 
the sphere, when calculated starting from: a) the true val- 
ues of fi\, up, n, adopted in the simulation (asterisks, solid 
line); b) the values of fj,\, fip, n, obtained after fitting the 
observations with a single star model (crosses, dotted line); 
and c) the values of fix, up, n and of the orbital parameters 
derived after fitting the observations with a 10 parameter 
model reproducing the standard motion plus the Keplerian 
circular motion induced on the observed star by the presence 
of the planet (diamonds, dashed line). 

The annual proper motion of 47 UMa is very large, com- 
pared to the magnitude of the astrometric signature pro- 
duced by the gravitational influence of its planetary com- 
panion: in ecliptic coordinates, — -0.283 arcsec/year, \ip 
= 0.152 arcsec/year, while a ~ 0.5 mas for an orbital in- 
clination of 45°. The top left panel shows the full range of 
motion of 47 UMa over five years; the looping motion is the 
parallactic eclipse, and the effect of the planet is essentially 
undetectable on this scale. In order to see the difference in 
the residuals, we zoom first on region A marked by a small 
square in the top left panel, then we identify in the upper 
right panel two 1-mas 2 regions (B and C) , representing each 
a lOOOx enlargement of a small fraction of the motion (bot- 
tom panels). The true and reconstructed motions are very 
close, consistent with the 10/xas single-measurement error 
of the observations used to derive the orbital parameters, 



while on this scale the orbit obtained with the single-star 
assumptions (dotted) is clearly very different from the ac- 
tual motion. 

Third, we derive the uncertainties in the fitted values of 
the important physical parameters of the two systems, and 
especially the companion masses. We derive the mass from 
the fitted orbital parameters via the mass function: 



i 

pa 



(5) 



{M s + M p ) 2 

where again M v and M 3 are the planetary and stellar masses 
in solar masses, it the parallax, P the orbital period in years, 
and a s the semi-major axis of the parent star, expressed in 
arcsec. 

Assuming M p <C M s , we obtain: 

1/3 



M v — I —^r 



»1 

r 3 p2 



(6) 



Fitting the simulated observations with a 10 parameter 
model we can derive estimates of tt and P directly, while 
a s can be obtained as a combination of the approximated 
four Thiele-Innes elements. Provided that one gets informa- 
tions about the mass of the primary, for example thanks to 
spectroscopy, it is possible to calculate M v . 

In our case we have, for simplicity, supposed to know 
M47UMa and Mi^vir to be exactly equal to M s . 

Tables [| and ^ show the approximated values of M v in 
the two cases, starting from the estimated values of P, n and 
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a 3 obtained during the fit, as function of orbital inclination. 
The Tables stop at i = 20° for 47 UMa and i = 10° for 
70 Vir, taking into consideration the u pper limits on their 
masses found by Perryman et al (1996) using HIPPARCOS 
data. The HIPPARCOS measurements, as a matter of fact, 
agreed with the stars to be single, as no evident signature 
greater than the nominal error (~1 mas) was revealed. This 
implies the following upper limits: M p ^. 7 Mj, for 47 UMa, 
Mp^i 38 Mj, for 70 Vir, which again means, considering the 
results obtained spectroscopically, choosing a minimum or- 
bital inclination as in the tables. In both cases, it is clearly 
evident the high accuracy with which the relevant orbital 
parameters and the planet's mass are recovered: thanks to 
the fact that S/N > 1, the fit to GAIA's simulated obser- 
vations of 47 UMa and 70 Vir is very satisfactory; the RMS 
errors between fitted and nominal values of the various pa- 
rameters are always well within 10 per cent. 

Such results confirm that GAIA could reveal itself a 
very powerful instrument to investigate many of the known 
candidate planetary systems, with the exception of very 
short-period systems, such as 51 Peg: provided the signal-to- 
noise ratio is sufficiently high and the period not too short, 
our simulations show that such systems will be easily de- 
tected, and their orbital parameters accurately determined 
by GAIA, at least if such systems are effectively simple, as 
the first spectroscopic measurements would suggest. In case 
the signatures produced cannot be interpreted as due to the 
presence of only one planet, this would introduce many in- 
terpretative complications in the signals observed: for plan- 
etary systems resembling our own, detection should be little 
affected and residuals analysis might reveal hints of the in- 
fluence of the outer planets (a for a Saturn-like companion 
is ~ 60 ^ias at 10 pc but its period is ~ 6 times the mission 
lifetime) . However, reliable orbital fitting would probably be 
restricted to the main component. On the other hand, mul- 
tiple orbital fitting to good accuracy should be possible for 
systems composed of giant planets of comparable mass, or- 
biting with periods within the interval ~0.5-5 yr, and yield- 
ing astrometric S/N > 10 , like it would be the case for the 
two o utermost planets of the recently discovered system v 
And flButler et al. 1999| ). 

The capability of detecting and measuring multiple 
planets with GAIA will become matter of future simulations 
and system studies. 



5 SUMMARY AND CONCLUSIONS 

In this work we have given the first quantitative evaluation 
of the detectability horizon of single extra-solar giant planets 
around single normal stars in the neighborhood of our solar 



estimates of the masses up to distances of order of 50 pc, 
value which doubles if we consider the range of periods in 
the vicinity of the mission lifetime; 

c) simulated observations of a selection of the actually 
known extra-solar planets, discovered by means of spectro- 
scopic measurements, provide a meaningful estimate of the 
uncertainty with which masses and orbital elements can be 
determined for the known star-planet systems and for a sub- 
stantial fraction of those that will be found within the con- 
text of such a global astrometry mission. Although prelimi- 
nary, our results indicate that these systems will be easy to 
discover and their orbital parameters will be accurately de- 
termined with GAIA, except very short-period systems such 
as 51 Peg. 



Hence, our results indicate that: 1) GAIA would moni- 
tor all of the hundreds of thousand F-G-K stars (i.e., whose 
masses are within a factor ~ 1.5 that of the Sun) up to a 
distance of ~ 200 pc from the Sun, in search for astrometric 
signatures due to the presence of giant planets (M ~ Mj) 
with orbital periods up to Jupiter's; 2) a significant frac- 
tion of the detected planets would have the main orbital pa- 
rameters (semi-major axis, period, eccentricity, inclination) 
measured to better than 30 per cent accuracy. 

Therefore, the GAIA survey would uniquely comple- 
ment the expectations from other ongoing and planned spec- 
troscopic and astrometric planet searches, both from ground 
and in space, thus helping with the creation of the fun- 
damental testing ground on which to measure the valid- 
ity of actual theoretical models of formation and evolution. 
GAIA's discovery potential might have significant impact on 
our knowledge of the distribution laws of the most relevant 
orbital parameters, and it would contribute to determine the 
frequency of planetary systems themselves in the solar neigh- 
borhood and, by extrapolation, in the whole Galaxy. A vast 
all-sky astrometric survey would help understand peculiar 
characteristics of these systems, e.g., whether giant planets 
lying in the outer regions are common: such planetary sce- 
narios may be worth further investigation, as, according to 
present theoretical models, this could indicate presence of 
low mass planets in the inner regions, possibly in the parent 
stars' habitable zones. The monitoring of hundreds of thou- 
sands stars directly implies the chance to investigate objects 
belonging to a wide range of spectral types, thus providing 
the important observational material (Boss 1998) to decide 



whether giant planets are more likely to form by means of 
gravitational instability in disks (once t hey are found more 



system i for the global astrometry mission GAIA. Complete 



simulat iono, comprohonaivo of oboorvationo of otar planet 



often around young s tars (Kuiper 1951; Cameron 1975; 



Bo- 



denheimer et al. 1980)), or by means of accretion of plan etes- 



system a aud successive statistical analysis of the simulated 
data, h ave yielded the following rooulto: 



imals (onc e they are found more often around o ld stars (Pol- 
lack 19841; iLissauer 19871 iPollack et al. 19961)). The high 



a) it will be possible to detect more than 50 per cent of all 
Jupiter-like planets (orbital period P = 11.8 years) orbit- 
ing solar-type stars within 100 pc; Jupiter-size planets, with 
shorter orbital periods, will be detectable up to 200 pc, with 
similar probabilities; 

b) for true Sun- Jupiter systems it will be possible to deter- 
mine the full set of orbital parameters and to derive accurate 



precision global astrometric measurements will estimate the 
inclination i of the orbital planes for the majority of the 
presently known planetary systems and for a large fraction 
of those that will be eventually discovered: it will then be 
possible to provide unambiguous mass estimations of such 
dark companions, reducing significantly the uncertainty on 
the mass range in the transition region from brown dwarfs 
to giant planets. 
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